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C°\ • Abstract. A cosmological model with a gravitational Lagrangian L g (R) oc R + AR n 

is set up to account for the presently observed re-acceleration of the universe. The 
evolution equation for the scale factor a of the universe is analyzed in detail for the 



two parameters n — 2 and n = 4/3, which were preferred by previous studies of the 



early universe. The initial conditions are specified at a red-shift parameter z ~ 0. 
The fit to the observable data fixes the free parameter A. The analysis shows that 
0^ the model with n = 2 agrees better with present data. Then, if we set w(q) = — 1 at 

z = 0, corresponding to the deceleration parameter q « —1/2, we find that at z « 0.5, 
w(q) has evolved to w ~ —0.72, corresponding to q ps 0. At z rs 1 we find u> ~ 
| corresponding to g sa 1/2. These results are compared with the flat Friedmann model 

with cold matter and Lambda-term (LCDM model) for the same initial conditions at 
z rs 0. The other choice of the model with n = 4/3 allows for big crunch. However 
this possibility is eliminated by the fit of A to the present data. 
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1. Introduction 

The discovery of a re-acceleration in the present expansion of the universe is interpreted 
as evidence for the existence of a cosmological A-term in the present equation of motion. 
At present, the most popular explanation is based on quintessence models which extend 
Einstein's action by scalar fields with different potentials. The additional parameters 
are chosen to comply with observable data. For a review see P3I2]- 

When studying such models it is important to match them with presently popular 
models of the early universe. Also here scalar fields have been used to modify Einstein's 
equations. But another important set of models makes use of curvature terms generated 
by fluctuating quantum fields as a Casimir effect, a mechanism discovered by Sakharov 

This gives rise to additional terms AL(R) which behave quite differently from the 
Einstein-Hilbert Lagrangian 

Starting from a scale invariant renormalizable action proportional to R 2 it is possible 
to generate for small R an Einstein action [6 plus terms of the form 

AL(R) = AR 2 + CR 2 In \R 2 /Rl\ (1) 

where A, C are parameters. The effect of such terms is studied with general methods 
developed an arbitrary function AL(R) in the Friedmann models [7 ■ 
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Such additional terms can, of course, be generated alternatively by an Einstein 
action coupled to a scalar field a with a suitably chosen potential [HE]- For instance, 
AL(R) = AR 2 is obtained from AL(R, a) = crR — a 2 /AA by extremization in a. 

There have been ample discussion of models in which AL(R) is chosen to fit 
observable data ^H]. Our choice ([TJ has the advantage having a physical origin and 
possessing only three parameters which can be fixed completely by present data at 
z = 0, thus allowing us to calculate uniquely the future and past evolution of the 
universe. 



2. Basic Equations 



As suggested by observations, we consider the flat cosmological Friedmann model with 
the metric 

ds 2 = dr 2 -a(r) 2 (dx 2 + dy 2 + dz 2 ). (2) 

If H denotes the presently observable Hubble constant (the subscript will always 
indicate the present epoch), the reduced curvature tensor p\ = HQ 2 R l k has the following 
matrix elements as a function of the reduced time 6 = H t, with the notation a = da/d6\ 

Pq = — 3d/ a, (3) 
pi = p = -6 (a/a + a 2 /a 2 ) . (4) 

The variation of Einstein's Lagrangian with an additional term AL(R) = f(R) gives 
Gt=^Tt + T?; G\ = P \-U\p. (5) 

H 2 
Here Zf corresponds only to cold matter in the present universe and 




W™) (Pi } (0) 



dp 



specifies the effective quintessence with the nontrivial dependence on curvature. 
We proceed as in Ref. ^T] by introducing the new variable 

V = (da) 2 (7) 



which allows to reduce the order of the equations. Then the i = k = component of 
Eq. © leads to 

( a dy\ a A df p 
\ y -2da -J fip) + ay ^ 



y + 



Pc -a 4 = n c a, (8) 



P* 

where p c is the a-dependent cold dark matter (CDM) energy density, p* = 3Hq/87iG is 
the critical density. Choosing the value of the scale factor a(9 ) equal to 1 at present, 
one has p c = po/a 3 and (p c / 'p*)a 4 = Q c a. 

In order to investigate the evolution of the cosmological model it is enough to obtain 
the solution of Eq. (jHj) with appropriate initial conditions. But for interpretation of the 
solution and for its comparison with observations it is necessary track for changes of 
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CDM energy density p c and quintessence energy density p v separately. For this purpose 
we use i = k = component of Einstein's equations (0) 

G° = P° -\p = ^r(p c + p v ). (9) 

Multiplying this equation by a 4 and using Eqs. (jHJ) and (J2J), we have 

y = (p c + p v )a i /p* (10) 

accounting for the evolution of cold matter energy density. From this we find for the 
quintessence energy density: 

Pl = (y-Q c a)/a\ (11) 
p* 

By solving Eq. (JHJ), we can find the evolution of p v . At known y, the Hubble parameter 
and the curvature components are obtained from 



We now use Eqs. © and (JIJ), rewritten as 

£ + # 2 = --, (13) 
6 a 

to derive an expression for the deceleration parameter of the universe 
da 

q = -— 14 

a 1 

in terms of the curvature components as follows: 

9 = ^ + 1- ( 15 ) 
For a comparison with observable data we relate the parameter q to quintessence 
and CDM energy densities. Let us suppose that the quintessence pressure and energy 
density are proportional to each: p v = wp v . Then the i = k = component of Einstein's 
equations becomes 

2 

: ) =HPc + Pv) (16) 



ay 

a) 

where K = 87rG and the equation for the scalar curvature gives 



6 



a /a x 2 



HPc + Pv~ 3p„) = N[p c + p v (l - 3w)]. (17) 



a \a 

From these equations we derive 

da 1 / 35 v w \ 

v = -^ = - 2 [ 1 + TTT v ) ] 5 v = PvIp- ( 18 ) 

The proportionality constant w is related to the deceleration parameter q by 
(2q-l\ {5 V + V 



The observed value of q at red-shift parameter 2 = a /a-l < 1 defines the values 
of y and y' Q in Eq. (jHj) in the present point a = 1. The data JTJ seem to imply that 
w < corresponding to the inflation stage. 
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3. Quintessence Model with / = Ap n 

Instead of our fluctuation-generated additional terms in a model with an additional 
term AL(R) = f(R) = Ap n oc R n has been considered in Ref. ^U]. Let us review this 
case briefly to see the difference with respect to our theory. The basic equation for dark 
energy from (jHJ) reads 




-2 



i -n /2 , fA sVy' 



n c a-y, (20) 



n(n — l)yy" H y' + (4 — 3n)n — 

2 a 

where 7 = (— 3) n_1 , y' = dy/da, where A is still arbitrary. Variation of the Lagrangian 
L(R) = f(R) gives, instead of Eq. ©, 

ft = KT t k , (21) 

where is energy-momentum tensor of usual matter. The authors in Ref. ^U] analyze 
these equation without taking into account of usual matter. In this approximation, the 
solution of our theory would be obtained by setting the tensor in Eq. (JUJ) equal to zero: 
T k = 0. For the choice f(R) = AR n , the corresponding equation is obtained by setting 
the brackets in Eq. (J2fl|) equal to zero. The solution is 



» = «»; k= . (4 "" 5)2 " ■ (22) 

(n — l)(2n — 1) 

By Eq. (|7J), this is corresponds to the following expression for the scale factor 

(n-l)(2n-l) . , 

a = 9 a ; a = -. 23 

2 — n 

This solutions has a constant deceleration parameter q. Further, to comply with the 
observable H Q , Ref. jTUj determined certain intervals of n with different ages of the 
universe. This allowed him to find the deceleration parameter q and parameter w (see 
Section as: 

g=(l/a-l); w = (2 - 3a)/3a. (24) 

Present data force us to reject such a pure power dependence of the scale factor in (J2BJ)- 
The latest data show that at present (z = 0), the universe is accelerating its expansion 
(w ~ —1) but at z = 1 there was no inflation (w ~ 0). Apart from that, the influence 
of cold matter cannot be neglected j!2j . 

Thus we need another ansatz for f(p) to account for present observations. From 
(f2T)j) we see that a simple power in the asymptotic solutions is absent for n = 2 and 
n = 4/3. These happen to be the same parameters which were of special importance 
in a previous theory of the early universe. At n = 4/3, one of us ^Hj has obtained for 
the first time a cosmological model without a singularity. That model passes through 
a regular minimum, had inflationary stage and tends asymptotically to the classical 
Friedmann solutions. In addition, the usefulness of an additional term R 2 in the models 
of the early universe has been pointed out before. Therefore we analyze the possibility 
of using of such powers for construction of models with variable parameters q and w. 
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4. Models with n = 2 and n = 4/3 

4.1. /(p) = -v 

This model was often used in the theory of the early universe. It describes the stage 
of fast oscillations of a{6) (the so called scalaron stage, which was introduced by A. 
Starobinsky in |14j). The damping of such oscillations was connected with creation of 
unstable particles and filling of the early universe by a hot plasma. 

Here we want to analyze this model in the opposite regime when the period of 
oscillations is commensurable with age of the universe 1/H . One can easily see the 
oscillations of the model by inserting the specified form of quintessence in (JHJ) and (JBJ) 

d 2 p ado , ^ . , , 

The scalar curvature performs oscillations near the value of p which corresponds to the 
model of cold dust matter in the Friedmann universe. 

More conveniently, we may use the i — k — component of the generalized 
gravitational equation (}20|) 



6A 



y"y - y' 2 /A - 2yy'/a] = -a 2 (y - Q c a). (26) 



Introducing the quadratic scale variable £ = a 2 , this becomes 



24A 



d 2 y 1 / dyV _ y_dy 
V de A\dd 2£d£ 



toJZ-V- (27) 



The initial conditions are specified at ao = 1, where the reduced Hubble parameter 
a'o/a is unity. This condition defines y = 1 at given £ = 1 [see Eq. (JJJ)]. 

At y Q = 1, the value of derivative (dy/da) is determined by observations in the 
following way. Using the observable estimation for Q v = 0.7 and Q c = 0.3, the value 
of 5 v q in Eq. (|18|) is approximately equal to 2. The observations for z << 1 show that 
w ~ — 1 and, using (|IB|L we find 

gb«i(l-2) = -l/2. (28) 

At £o = 1, we obtain from (|18|) 

(tZy/de) w 3/2. (29) 

The only open parameter left is A in Eq. ()27j) . It specifies the amount of dark energy 
in the universe. 

The solution y(£) in the interval 0.01 < £ < 1 is presented in Fig. 1. The age of 
the universe T is defined by the value of the Hubble parameter 

~ = KZ) = y/m/Z, 8 = TH = f-¥=, 6^0.(30) 



Note that Eq. (|2*7j) yields for £ — > the regular behavior y — > const. Thus the integral 
at e — > has no singularity. 
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Figure 1. For the models with n = 2 (left graph) and with n = 4/3 (right 
graph) the following parameters subject to £ are shown. 



The age T depends on choice of the parameter A in ()27|) . The satisfactory value 
T pa 0.83/iJo ~ H-4 Gyr is reached at A ~ 1. For a further analysis of the model and 
comparison with observable data see Section 03 

12. f(p) = -Afl* 

The special feature of the model with n = 4/3 consists in absence of the scale factor a 
in explicit form in the Einstein's equations if matter is neglected jTSj. This allows one 
to find the general solution of Eq. (|20J) with given n. One can show that de Sitter's 
solution arises in the limit a ^> 1. The curvature of this limiting solution is determined 
by the parameter A. 

The evolution equation for this case is from (J2U|) : 



y- 



d 2 y 




+ 



lydy 



9 1 



'2dy\ 2/3 

Add 



^c\r - y 



(31) 



<% 2 8\d£) ' 2£d£ 16 
By analogy with previous subsection, we obtain for this model (see Fig. 1) an age 
of the universe T w 0.79/ H pa 10.8 Gyr with a parameter A = 0.33. 

Leaving the analysis of the solution to Section let us consider here only the 
following issue. For £ > 1, the influence of cold matter becomes negligible. Then 
Eq. (|31[) has a first integral. It can be easily seen from (|2Uj) for n = 4/3: 

2 / , \ 2/3 , 





2/ j 



(3)- 1/3 /A 



(32) 
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The first integral is 




) 



4/3 



-D(y + C^y) 



(33) 



where C denotes an integration constant. For C > 0, this equation has the asymptotical 
solution y pa 7a 4 , 7 = 2(D/3) 3 , which coincides with the de Sitter solution 



For C < 0, the solution possess a scale parameter a where dy/da = 0. Beyond that, 
y starts decreasing. The solution reaches a point y — > 0. By definition, 7/ can never be 
be negative. Analytic continuation of the solution beyond the specified point leads to 
an increase of y at decreasing of a. As a approaches zero we obtain a new singularity 
[7j. The zero of y corresponds to a maximum of the scale factor a with subsequent 
contraction of the universe ending in a singularity. 

This is an example for a big crunch (see £Q and references therein). It is important 
to note that regime of a crunch is realized here for a flat cosmological model. 

Note that a contraction of cosmological model to a point cannot really be described 
by a homogenous model. In the limiting regime, the universe becomes certainly 
inhomogeneous and anisotropic because of fast growth of perturbations. Using of the 
solution of Eq. (|3*3*|) with parameters D and C corresponding to observable data for 
z < 1 guarantees an absence of a crunch in future (C > 0). 

5. Discussion 

After fitting the model's parameters to the present observable data - acceleration of 
the universe, the Hubble parameter at red-shift parameter z — 0, and the age of the 
universe - there are no free parameters in the model. Its predictions for large z can be 
compared with observations. 

As can be seen in Fig. 1 (in case of n = 2), the acceleration of the expansion 
changes at z ~ 0.46 to a deceleration (w ~ —0.72, q pa 0). Near z = 1 the variant of 
dust-like dark energy (w = 0) is realized. It corresponds to latest observational data [I]. 

A special feature of the model is that for a — > 0, the variable y = (da) 2 tends 
to a constant, which is equivalent to an evolution of the universe filled by hot matter 
(p = e/3) in the Friedmann cosmology. When the solution is oscillating, the inflation at 
z pa cannot be eternal although the period of oscillations is comparable with the age 
of universe. 

In case of n = 4/3, the crossover from a decelerated expansion to an accelerated 
one takes place at z pa 0.30 when w ~ —0.60. 

It is interesting to compare our model with the simplest LCDM model (for review 
see PP). Using the notation of Section (J2J), we obtain 





(34) 



y = 1U V2 + ^U 2 



(35) 



with the same values f2 c ~ 0.3 and f2. 



'V 



0.7 
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Our results show the following: the age of the universe is near 1/Hq ~ 13.7 Gyr. 
This age is larger than the age from the pure p 2 -model (~ 11.4 Gyr) and from the 
p 4 / 3 -model (~ 10.8 Gyr). The parameter w of the LCDM model calculating with use of 
()19|) has the value —1. We may use Eq. ()18|) to determine the deceleration parameter 
q and compare it with the p 2 and p 4 / 3 models. At z — 1, the deceleration parameter is 
q ~ 0.2, whereas the p 2 and p 4 / 3 models have q ~ 0.5. This value is closer to observable 
data. At z ~ 0.5, the parameter q in the p 2 -model is close to zero and in the LCDM 
model q ph —0.1. Apparently, the p 2 -model corresponds better to observations. The 
LCDM model is mathematically very simple. But it leaves the value of the A-term an 
unexplained fundamental constant. For dynamical models of the A-term (for example, 
the p 2 -model) this value evolves from a large Planck value in the early universe to small 
value at present. 
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